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ABSTRACT 

We present new limits on cosmological parameters from the statistics of gravitational 
lensing, based on the recently revised knowledge of the luminosity function and internal 
dynamics of E/SO galaxies that are essential in lensing high-redshift QSOs. We find 
that the lens models using updated Schechter parameters for such galaxies, derived 
from the recent redshift surveys combined with morphological classification, are found 
to give smaller lensing probabilities than earlier calculated. Inconsistent adoption of 
these parameters from a mixture of various galaxy surveys gives rise to systematic 
biases in the results. We also show that less compact dwarf-type galaxies which 
largely dominate the faint part of the Schechter-form luminosity function contribute 
little to lensing probabilities, so that earlier lens models overestimate incidents of 
small separation lenses. Applications of the lens models to the existing lens surveys 
indicate that reproduction of both the lensing probability of optical sources and the 
image separations of optical and radio lenses is significantly improved in the revised 
lens models. The likelihood analyses allow us to conclude that a flat universe with 
Uo = 0.3i/gand flo + V = 1 is most preferable, and a matter-dominated flat universe 
with Ao = 0 is ruled out at 98 % confidence level. These new limits are unaffected by 
inclusion of uncertainties in the lens properties. 


Subject headings: cosmology: observations — gravitational lensing — galaxies: general 


1. INTRODUCTION 

Recent cosmological observations have increased various lines of evidence that the matter 
density of the universe, Uo, falls short of the closure density required to halt its expansion (Ostriker 
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& Steinhardt 1995; Yoshii & Peterson 1995; Krauss 1997; Bahcall & Fan 1998). These include 
the anisotropies of the cosmic microwave background, large-scale structure of galaxy distribution, 
number evolution of clusters of galaxies, and large amount of gas in clusters combined with the 
constraints on light-element nucleosynthesis. Then, in order to reconcile with the flat geometry 
of a universe favored by inflationary models, a cosmological constant A, or Ao = Ac 2 /377q where 
Hq is the Hubble constant, has been invoked to achieve the critical energy density, Ho + Ao = 1. 
In fact, the existence of a non-vanishing cosmological constant points to a concordance of various 
observations and cosmological models with a combination of Ho ~ 0.3 and Ao ~ 0.7 (Ostriker 
& Steinhardt 1995). The case for a flat universe with a non-zero Ao appears to be more likely 
than an open universe, even when the recent revisions of the mean age of globular clusters and 
the current estimate of Hq are incorporated (Krauss 1997; Pont et al. 1998). Most recently, the 
distance determinations of Type la supernovae at cosmological distances (Perlmutter et al. 1998; 
Garnavich et al. 1998) revealed further evidence for an accelerative expansion of a universe. 

On the contrary, tight limits on a cosmological constant have been put from the statistics 
of gravitational lensing, because the number of multiply imaged QSOs found in lens surveys is 
sensitive to the value of Ao (Turner 1990; Carroll, Press & Turner 1992). Systematic surveys 
to search for lensed QSOs, such as the Hubble Space Telescope (HST) Snapshot Lens Survey 
(Maoz et al. 1993), have identified only a few lenses among hundreds of QSOs, and this small 
lensing probability has been the central argument against a non-vanishing cosmological constant 
(Maoz & Rix 1993, hereafter referred to as MR). In particular, Kochanek (1996) put severe limits 
Ao < 0.66 at 95 % confidence in flat cosmologies, using all of existing optical lens surveys available 
then, together with the data of radio lenses to reproduce their image-separation distributions. 
These constraints on Ao challenge the proposed cosmic concordance required from other various 
observations. 

Modeling the statistics of gravitational lensing however contains uncertainties associated with 
the dynamical structure and number density of galaxies that work as a lens. Among various galaxy 
types, elliptical (E) and lenticular (SO) galaxies play an essential role in the models, because 
spiral and irregular galaxies are known to make negligible contributions to lensing statistics (MR; 
Kochanek 1996). Thus, attempts have been made to tighten the lens properties from spectroscopic 
and photometric data of E/SO galaxies (Fukugita & Turner 1991, hereafter FT; Breimer & Sanders 
1993; Franx 1993; Kochanek 1994). The internal dynamics of such galaxies is deduced from the 
velocity dispersions and light profiles, whereas the luminosity function (LF) of E/SO galaxies 
provides the number density of galaxies responsible for lensing background QSOs. Thus, the 
resultant constraints on Ao inevitably depend on how well all of these lens properties are modeled. 

Here we present new calculations of the lensing statistics in view of the recently updated 
knowledge of the E/SO-type LFs and their internal dynamics. We show that the models using 
the updated Schechter parameters for E/SO-type LFs, which have been derived from the recent 
redshift surveys combined with morphological classification, significantly revise the theoretical 
predictions of both the lensing probabilities and image-separation distributions. We also point 
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out that earlier models considerably overestimate the lensing effects of less compact dwarf-type 
galaxies, which dominate the faint part of the Schechter-form LF. Our calculations using realistic 
internal structure and dynamics of such galaxies in §3 are found to diminish their contribution 
to lensing statistics. Our revised calculations applying to the data of the HST Snapshot Lens 
Survey were reported in Chiba & Yoshii (1997), showing that a low-density, flat universe is not 
statistically inconsistent with observations. This article is devoted to the further properties of 
our models for the lensing statistics in more detail. In particular, in order to pursue the most 
plausible set of cosmological parameters (Ho,Ao), we use all °f the optical lens surveys including 
the HST Snapshot Lens Survey. The data of radio lenses are also taken into account to tighten the 
results of the statistics (Kochanek 1996), by estimating the relative likelihood of their observed 
image-separation distributions. In §2 we present the models for gravitational lensing and set all 
of ingredients required in the lensing statistics. The effects of the E/SO-type LFs on the lensing 
probability are discussed in §3 to demonstrate how the updated knowledge of the LFs modifies 
the results in a significant manner. Then in §4, we apply the current lens models to the existing 
lens data and derive the new limits on cosmological parameters. Implications of our results are 
discussed and conclusions are drawn in 85. 


2. GRAVITATIONAL LENSING OF QUASARS 


2.1. Lens model and lensing probability 


In the presence of a lens along the light-ray path from a source to an observer, the ray is 
deflected by an amount that depends on the gravitational field of a lens. We assume that the 
density profile of a lens is represented by that of an isothermal sphere having a finite core, 


p(r) 


2vrG(r 2 + r 2 core ) ’ 


( 1 ) 


where a is the one-dimensional velocity dispersion and r core is the core radius. The impact 
parameter of a deflected light in the lens plane, b, is related to that of an undeflected light, l, 
through the equation (Hinshaw & Krauss 1987), 


b 

where r e is the critical impact parameter for multiple imaging in the case of a singular isothermal 
sphere (SIS) as a lens. This is defined as r e = 4ir(cr/c) 2 DolDls/D os, where Dol is the 
angular-diameter distance between the observer and the lens, Dls between the lens and the 
source, and Dos between the observer and the source. We use the filled beam formulae to define 
these distances. 


Equation (2) admits three solutions bi(l) with i = 1,2,3 if l is less than a critical impact 
parameter Iq [eq.(A2)-(A4)], under the condition that the ratio between r core and r e , (5 = r core /r e , 
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is smaller than 1/2. The image 1 lies on the same side from the lens as the source, whereas the 
other two images (images 2 and 3) lie on the opposite side of the lens from the source, where 
the innermost image 3 is faintest among three images. If either l > Iq or (3 > 1/2 is fulfilled, 
there is only one image. We summarize the characteristic properties of these lensed images in 
the Appendix (see also Hinshaw & Krauss 1987). Each image is then magnified by a factor 
Ai = \bidbi/ldl\ compared to an unlensed image, and the total magnification factor A is defined as 
A = A\ + A2 + A 3 . As demonstrated in the Appendix, the presence of a core leads to larger A 
than the SIS case, while the cross section for lensing is reduced, ttIq < -nr 2 . The image separation 
angle Ad for these multiple images is well characterized by that between outer two images (images 
1 and 2 ) for l = 0 , namely, A 9 = 2(1 — (3) l / 2 r e /D ol- 

Statistics of gravitational lensing depend on the number density distribution of lens galaxies. 
This is deduced from the LF of galaxies 1 j> g (L g ) with luminosity L g . To compare with the earlier 
standard models of lensing statistics, we adopt the assumption that the comoving number density 
of galaxies is independent of lens redshift zl- The results of the Canada-France Redshift survey 
by Lilly et al. (1995) support this assumption (see also Totani & Yoshii 1998). 

Then, the probability that a QSO with redshift zs and luminosity Lq is lensed is evaluated 
from the fraction of the QSOs that are amplified to this luminosity by lensing, among the total 
number of QSOs. Using the LF of QSOs defined as 4>q(Lq), this is expressed as, 


p(Lq,z s ) = 


f ZS , /, 

/ dz L (l + z L f -— 

/0 dZL . 


dLg(f)g {Lg) 


r ‘ 0iL ’' ZL \i2nis ^sIAAA 

0 (Pq{Lq,zs) 


( 3 ) 


where S denotes the selection function which depends on the image separation and magnitude 
difference between the primary and secondary images, and on the seeing in a lens survey (e.g. 
Maoz &; Rix 1993). To perform the integration over l, it is convenient to avoid the divergence of 
image magnification when a QSO lies at l = 0 and Iq. For this purpose, we integrate it over the 
position of image 2, 62 , instead of the source position l, using the Jacobian | ^ 2 1 which relates these 
two positions (Blandford & Narayan 1986). This yields / dl2irl —> / <^ 6227 T &2 /1 ^2 1 over the range 
of b- < b 2 < b + , where b + and b _ are the critical positions of image 2 relevant to l = 0 and Iq, 
respectively. 

In order to relate luminosities of galaxies L g to their structural parameters (a, r core ), we first 
assume the scaling relation between L g and a, L g = L* g {a/a*) 1 , where er* is the one-dimensional 
velocity dispersion at L g = L*. This is in analogy to the Faber-Jackson relation for E and SO 
galaxies, and we adopt 7 = 4 for these galaxies that effectively work as a lens. As discussed 
in Kochanek (1996) and Chiba & Yoshii (1997), we will not introduce the correction to a* for 
the dark-matter velocity dispersion by a factor of (3/2) 1 / 2 , because detailed dynamical modeling 
of early-type galaxies by Breimer & Sanders (1993), Franx (1993), and Kochanek (1994) has 
invalidated this correction. If this correction is employed, the models will overestimate the lensing 
probability by a factor of 2.25. Second, in contrast to the simple treatment of r core = const in 
Chiba & Yoshii (1997), we take into account the observed increase of core radii with increasing 
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luminosities of E/SO galaxies, by the relation L g = L*(r core /r core where 77 = 1.2 (FT; 
Kochanek 1996; Faber et al. 1997). From these scaling relations, the parameter (3 = r core /r e is 
expressed as /3 = {r* ore /rl)(L g /L*) 11 ^ 1 / 2 , thereby increasing with increasing L g . 

In conventional lens models, the LF of lens galaxies 4> g is assumed to hold the Schechter form, 

ML 9 )dL g = r 9 (L 9 /L* g ) a e-( L ^dL g /L* g , (4) 

where </*, a , and L* are the normalization, the index of faint-end slope, and the characteristic 
luminosity, respectively. In this case, the lensing probability p is proportional to a non-dimensional 
factor F*, defined as 

F*^167r 3 (^) 4 (-|-) 3 ^r(a + 4/ 7 + l) , (5) 

(Turner, Ostriker & Gott 1984; FT). 


2.2. LF of QSOs and magnification bias 

Equation (3) indicates that the lensing probability depends on the shape of the QSO LF, 
4>q{L q ). As a standard case, we use the ih-band LF adopted by Wallington & Narayan (1993) 
(hereafter WN) on the basis of the LF determination for zs < 3 (Boyle et al. 1988) and that 
for zs > 3 (Schmidt et al. 1992; Warren et al. 1992; Irwin et al. 1992). Written as a function 
of absolute magnitude, Mq, the QSO LF = \dLq/dMffiypq is given in the smoothed two 
power-law form 

$q(Mq,Zs) oc 10 o. 4( Mq _ Mz )( j3i+1 ) + 10 o.4(Mq-M*)(/3 2 +1) ’ ^ 

where M z is the magnitude of the break in the LF. For zs < 3, this is given as 

M z = -2.5 Ll log(l + zs) + M 0 , (7) 

with ki, = 3.15 and M 0 = —22.421? mag for h = 0.5 (h = 77o/100 km s _1 Mpc^ 1 ), = 1, and 

Ao = 0. For zs > 3, the evolution of the break magnitude M z is corrected by — 0.54(zs ~ 3) 
to accord with the high-redshift QSO surveys (Wallington & Narayan 1993). The parameters 
Pi = —1.44 and = —3.79 denote the faint and bright end slopes, respectively. 

The effect of the so-called magnification bias such that lensed QSOs are over-represented 
in magnitude-limited sample due to image amplification is taken into account through the ratio 
c))q{Lq / 'A, zs) / 'A/ (/)q(Lq, zs) in the integrand of eq.(3). In contrast to the SIS case, we note that 
in the presence of a finite core, the integration over l (or over 62 ) that includes this ratio cannot 
be taken outside of other integrations over z^ and L g due to the dependence of A on P(zi,L g ) 
(see also Harnana et al. 1997 for more detail). To demonstrate the amount of the biasing effect, 
designated as B(Lq,zs), we use the following definition 

V 

p(A -> 1 ) ' 


B(Lq,z s ) 


( 8 ) 
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In the SIS case, this is reduced to B = f£° dA(8/A 3 )(j)Q(LQ/A, zs)/A/(/)q(Lq, zs ). 

We calculate the bias B(Mq,zs ) as a function of Mq in the B band, using a set of fiducial 
lens parameters 7 = 4, 7] = 1.2, and a* = 220 km s^ 1 . Figure 1(a) compares the SIS case 
(solid lines) with the case having a finite core r* ore = 0.1/U 1 kpc for a = — 1 (dotted lines) 
and a = 0 (dashed lines). It is evident that the bias is made large in the presence of a core, 
so that the effect of smaller cross section for multiple imaging than in the SIS case is reduced 
(Kochanek 1996; Hamana et al. 1997). This is attributed to the larger image amplification 
A in the presence of a core. It is also noted that the shallower faint-end slopes a in the LFs 
of galaxies yield the larger bias, which thereby reduces the effect of less numbers of galaxies 
in such LFs on the lensing probability. This arises from the fact that the LFs having the 
shallower a weight toward larger L g galaxies, thus larger 0, while the magnification bias increases 
with increasing (3 (Hamana et al. 1997). Figure 1(b) shows the effects of changing QSO 
LFs on the bias. In comparison with WN’s QSO LF, we also consider the broken power-law 
QSO LFs adopted by MR (0\ = — 1.2,02 = —3.6, = 3.5, Mo = —21.75 mag), and FT 

(Pi = —1.7,02 = —3.15,break in apparent magnitude = 19.15 mag). When compared to WN’s 
LF, FT’s LF yields the significant reduction of the bias at the bright end due to its much shallower 
slope 02, whereas MR’s LF yields roughly the same bias although it is slightly smaller at high 
redshifts zg. This suggests that the lensing probability based on either FT’s or MR’s LF will be 
smaller than that based on WN’s LF, as was also noted in Kochanek (1996). 


3. E/SO-TYPE LUMINOSITY FUNCTIONS AND LENSING PROBABILITIES 

3.1. Effects of updated Schechter parameters 

Many of earlier lensing statistics use the Schechter-form LF of all-type galaxies in Efstathiou 
et al. (1988) scaled by the Postman & Geller (1984) fraction (31%) of E/SO galaxies, because these 
early-type galaxies are essential in lensing statistics (FT; MR). The corresponding E/SO-type LF 
is hereafter referred to as EEP’s LF and its Schechter parameters (M* ,a,(f*) are tabulated in 
Table 1, where M* is the characteristic magnitude derived from L*. Although this procedure of 
counting the number density of E/SO galaxies has widely been adopted for lensing statistics, the 
different sampling and calibration between these observational studies prevent us from knowing 
the detailed shape of the E/SO-type LF in the sample of Efstathiou et al. (1988). Specifically, the 
faint-end slope a and characteristic magnitude M* for only E/SO galaxies are not necessarily the 
same as those for all-type galaxies. Thus, it is not surprising that the lensing statistics based on 
EEP’s LF holds some unavoidable systematics. 

Comparing with the above indirect method, the type-dependent LF for only E/SO galaxies 
has been obtained on the basis of the recent redshift surveys of field galaxies combined with 
morphological classification. We consider the two LFs derived from the Stromlo-APM survey 
(Loveday et al. 1992, LPEM) and the CfA survey (Marzke et al. 1994, MGHC), in addition to 
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EEP’s LF. The Schechter parameters for only E/SO galaxies are taken from these papers and 
tabulated in Table 1. Comparing with EEP’s LF, LPEM’s LF is characterized by the shallow slope 
a at the faint end and the small normalization (/>*, whereas for MGHC’s LF, M* is about 1 mag 
fainter and cf>* is about twice larger than the respective values of EEP’s LF. These revised values 
of the Schechter parameters will yield the different values of a* and F* from those based on EEP. 

We note that the high normalization </>* of MGHC’s LF may be explained in terms of local 
density inhonrogeneities present in the CfA survey volume, while the Stromlo-APM survey samples 
a volume ~ 30 times larger than that of the CfA survey and so provides the reliable space density 
of galaxies. However, the reason for the faint value of M* in MGHC’s LF is not well understood. 
There are some discussion that this is due to the fainter Zwicky magnitudes used in the CfA 
survey than the bj magnitudes in the Stromlo-APM survey, by an amount of 0.2 — 0.5 mag (see 
e.g. Huchra 1976; Auman et al. 1986; Bothun & Cornell 1990, about detailed discussions of the 
Zwicky magnitudes). If we transform M* in MGHC by adding a constant offset, then the resultant 
value of M* is shifted toward the values in EEP and LPEM. To see how this offset changes the 
lensing probability, we also consider the case of M* —» M* + 0.5 mag in MGHC, hereafter referred 
to as MGHC2 (Table 1). We note that there are also some suggestions that the Zwicky magnitudes 
depend on surface brightness of galaxies, in a manner that the luminosities of compact galaxies 
such as E’s tend to be overestimated (Auman et al. 1986; de Vaucouleurs et al. 1991). If this 
is the case, M* for E galaxies measured from the Zwicky magnitudes may be brighter than from 
other magnitude systems. As discussed in EEP, the transforms between the magnitude systems 
also contain color terms, complex isophotal corrections, etc., so that the procedure of transforming 
M* from one system into another by adding a constant offset is a gross oversimplification. Thus, 
in view of these uncertainties, the results based on MGHC2 need caution. 

Most recently, the Las Campanas Redshift Survey (LCRS) provided the expanded galaxy 
sample in the largest survey volume (Lin et al. 1996). The LF of all galaxy sample (after 
transformation from the Kron-Cousins R band system used in LCRS to the bj band) is in good 
agreement with that from the Stromlo-APM survey. Although the type-dependent LFs are not 
yet available, the LFs based on the [O II] A3727 emission feature (Lin et al. 1996) or on the 
more detailed spectral classification scheme (Bromley et al. 1997) show that the faint-end slopes 
a become progressively shallower from late to early-type galaxies. This dependence of a on 
galaxy type is in good agreement with LPEM’s result. Thus, an argument that the shallow a in 
LPEM’s E is due to failure of identifying faint ellipticals (MGHC) may not be supported. Here 
for comparison, we also take into account LCRS’s LF based on the galaxies that show no [O II] 
A3727 emission (Table 1). We shift LCRS’s M* by 1.1 mag to match the mean rest-frame color 
< bj — R> o= 1.1 of LCRS galaxies (Lin et al. 1996). It is worthwhile to note that the derived M* 
is slightly fainter than those in other LFs. However, the results also need great caution because of 
ambiguous morphological classification. 

In order to derive the characteristic velocity dispersion a* from the observed M*, we use the 



Faber-Jackson relation determined by de Vaucouleurs & Olson (1982) for early-type galaxies, 

— M* + 5 log h = 19.37 + 10(log a E — 2.3) for E 
— M* + 5 log h = 19.75 + 10(log a E0 — 2.3) for SO . (9) 

Then, for the fractional values (/e, fso) of (E, SO) galaxies in the number density, respectively, we 
evaluate a* for E/SO galaxies by means of a* 4 = {/e^*e 4 + f so^sq 4 ) / {fE + fso)- This is because 
the lensing probability is proportional to the fourth power of a* in the definition of F* (eq.5). 
We adopt the Postman & Geller (1984) fraction (E:12%, S0:19%) to calculate a* for the LFs of 
EEP, LPEM, and LCRS, while for the LFs of MGHC and MGHC2 we use Marzke et al. (1994)’s 
Table 1. The values of cr* and F* derived from these new LFs appear to be systematically smaller 
than those from EEP’s LF (see Table 1) due to the above mentioned difference in the Schechter 
parameters. Therefore, the total lensing probability based on the updated Schechter parameters 
for E/SO galaxies will be reduced. 

To highlight these effects of updating the Schechter parameters, we calculate the lensing 
probability of a QSO with Mq = —25.5 + 5 log h mag and zg = 2, for the standard lens parameters 
of 7 = 4, r* ore = 0.1 hr 1 kpc, and r/ = 1.2. Figure 2 shows the image-separation distributions 
p(Ad) for the cases of EEP (solid line), LPEM (dotted line), MGHC (short-dashed line), MGHC2 
(long-dashed line), and LCRS (thick solid line). The image-separation distribution for the case of 
EEP differs significantly from those for LPEM and MGHC. LPEM’s LF having shallower faint-end 
slope a and smaller normalization </>* shifts the distribution to somewhat large image separations 
with smaller overall amplitude. MGHC’s LF having fainter M* considerably suppresses the 
distribution at larger separation and makes it peaked at smaller separation. In either case, the 
total lensing probability is decreased when compared to the result from EEP’s LF. This clearly 
indicates that the predicted lensing probability is sensitive to the adopted Schechter parameters 
for the LFs of E/SO galaxies. The figure also suggests that the case of MGHC2 produces a similar 
p(A0) to that of EEP, whereas the case of LCRS leads to the skewed p(A9) towards larger image 
separations with smaller overall amplitude. Considering great ambiguities in the E/SO-type LFs 
inferred from MGHC2 and LCRS as discussed above, we will not consider these cases in what 
follows. 

Updated LFs of LPEM and MGHC hold separate combinations of Schechter parameters 
and yield different lensing probabilities. However, Kochanek (1996) contrived the LF in a mixed 
manner by adopting M* from the all-type LF by EEP, a from the all-type LF by LPEM and 
MGHC, and 4>* from the all-type LF by LPEM after scaled with MGHC’s fraction of E/SO 
galaxies. In his lens model, the characteristic velocity dispersion a* for E/SO was taken from the 
work of Kochanek (1994) that showed a* = 225 km s -1 using EEP’s M*. To compare with other 
cases, we also consider the LF adopted by Kochanek (1996), hereafter referred to as K96’s LF, in 
the statistics discussed in §4, and the associated Schechter and lens parameters are given in the 
last line of Table 1. We remark that in each galaxy survey, the Schechter parameters have been 
determined in a highly correlated manner, so that Kochanek’s method of combining the results 
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from various surveys lacks consistency and thus causes artificial systematics. In particular, while 
the number density (j>* of E/SO galaxies was averaged between LPEM and MGHC, adoption of 
the bright M* = —19.9 mag and steep a = —1, irrespective of the reported fainter M* in LPEM 
and the shallower a in MGHC for E/SO galaxies than each other, leads to a bias in favor of large 
lensing probabilities. The best way of avoiding these artificial systematics is to use the Schechter 
parameters as reported in the original references. 


3.2. Effects of deviation from the Schechter form 

Earlier models for lensing statistics adopted the Schechter form for the E/SO-type LF. 
Specifically, the index of faint-end slope a is usually ~ —1.0 as in EEP, thereby implying the 
existence of numerous faint E/SO galaxies that work as a lens. However, the faint part of the 
E/SO-type LF is still uncertain because the limited resolution hinders reliable type classification 
for faint galaxies. Moreover, distance determinations for field galaxies from their observed redshifts 
may be uncertain to some extent because of the peculiar galactic motions in the local field. On the 
other hand, the problem of accurate distance determinations is avoided by using nearby clusters 
of galaxies which offer a homogeneous sample of member galaxies at the same distance, observed 
down to faint magnitudes. Sandage, Binggeli & Tammann (1985) and Ferguson & Sandage (1991) 
have indicated that while the faint part of the LFs for all galaxies is dominated by diffuse and less 
luminous dwarf-type galaxies, the LFs for E/SO galaxies alone also show a sharp decline towards 
faint magnitudes, similarly towards bright magnitudes. This bounded form of LFs at both bright 
and faint ends was first recognized by Hubble (1936). For field E/SO galaxies, there is also an 
indication of such a decline towards faint magnitudes (Binggeli, Sandage & Tammann 1988; Driver 
et al. 1994; Driver et al. 1995). 

It has also become clear that early-type dwarf galaxies below M = —17 ~ —181? mag which 
populate at the faint end of the LF show the different internal structure and dynamics from 
those of bright ellipticals (e.g. Ferguson & Binggeli 1994). The radial brightness profile of such 
dwarf galaxies is well represented by the exponential form I(r) oc exp(— r/r^), where r^ is the 
scale length. This is reminiscent of the profile of spiral disks and is in contrast with the more 
centrally concentrated r 1 / 4 law in bright ellipticals. This brightness profile suggests a large core 
radius r core — 1 kpc (e.g. Kormendy 1988; Ferguson & Binggeli 1994), so that the parameter 
/3 = r core /r e introduced in the current lens model may be larger than 1/2 in all lens redshifts. 
Also, the luminosity-velocity relation follows L oc cr 2 ' 5 (Held et al. 1992) that deviates from the 
L oc a 4 law of the bright E/SO galaxies. Thus, the faint galaxies have smaller internal velocity 
dispersions than that expected from the L oc a 1 law. These diffuse properties of faint early-type 
galaxies are well predicted from the models of formation of dwarf galaxies invoking galactic winds 
and subsequent adiabatic expansion (Dekel & Silk 1986; Yoshii & Arimoto 1987). Consequently, 
because these faint early-type galaxies have large r core and small cr, their contribution to the 
lensing statistics is negligible. 
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To demonstrate these properties of dwarf galaxies in lensing calculations, we construct the 
model LFs of early-type galaxies which consist of both E/SO and dE/dSO types. The parameters 
for the LFs used are tabulated in Table 2 and the shape of the LFs is graphically displayed in 
Fig.3. The Schechter parameters for E/SO-type galaxies roughly accord with those derived from 
existing galaxy surveys, except that the LF at the faint magnitude is reduced by multiplying 
a factor exp(—io°- 4 ( M_A ^ cu )) with M cut = —17 mag in order to match the observed shape of 
E/SO-type LFs (Driver et al. 1995). The Schechter parameters for dE/dSO-type galaxies are taken 
from the work of Phillips and Driver (1995), by assuming M* (dE/dSO) = M*(E/S0 ) + 3 mag, 
<j>* g {dE/dS 0) = cj)*(E/S 0), and a very steep slope a = —1.5. As is evident from Fig. 3, the total 
LF (thin solid line) derived from the E/SO-type (dashed line) and dE/dSO-type LFs (dotted line) 
resembles a single Schechter form. The parameters by which this total LF is fitted to a single 
Schechter function (thick solid line) are also given in Table 2. 

Using these LFs, we calculate the lensing probability of a QSO with Mq = —25.5 + 51og/i 
mag and zs = 2. For E/SO galaxies, the lens parameters are taken as 7 = 4, r* ore = 0.1/r -1 kpc, 
and rj = 1.2, whereas for dE/dSO galaxies, we assume 7 = 2.5, r core = 0.2 h ^ 1 kpc ~ constant 
( 7 / = 0). In order to evaluate er* from M* for E/S0 galaxies, we use eq.(9) and the Postman & 
Geller (1984) fraction for these galaxy types. For dE/dSO galaxies, we use the data of IC794 
(Bender & Nieto 1990; Held et al. 1992) to obtain the following approximate relation that obeys 
the scaling L g /L* = (a/a*) 2 ' 5 , 

- M* + 5 log h = 16.00 + 6.25(log cr* - 1.72) for dE/dSO . (10) 

Figure 4(a) shows the lensing probability at zl = 0.5 as a function of galaxy magnitudes M 
in the B band. It follows that dE/dSO galaxies contribute very little to the lensing probability; if 
we adopt a more realistic core radius of r core ~ 1 kpc (e.g. Kormendy 1988; Ferguson & Binggeli 
1994), the lensing probability is much more reduced than presented with r core — 0.2 kpc. Figure 
4(b) shows the image separation distribution. It is evident from these figures that use of a simple 
Schechter form for E/SO-type LF, as has been adopted in previous lensing models, considerably 
overestimates the lensing probability at the faint magnitudes and therefore at the small image 
separations. 

In order to consider this trivial effect of dwarf-type galaxies on lensing probabilities at 
the faint part of the LFs of EEP, LPEM, and MGHC, we introduce a magnitude cutoff of 
exp(— iQ 0A ( M ~ MCVt )) in these LFs and examine the effect of changing M cut over its possible range. 
Figure 5(a) shows the image-separation distributions p(A9) for M cut = +00 (no cutoff), —16, and 
— 17 mag, and Fig. 5(b) shows the ratios of the total lensing probability p with cutoff relative to 
that without cutoff as a function of M cut . It is clear that the suppression of faint, less-massive 
galaxies from the originally Schechter-form LF reduces the incidents of small-separation lensing. 
We also find that the result from LPEM is very insensitive to M cut , because the shallow faint-end 
slope a of LPEM’s LF implies the small number of faint E/S0 galaxies irrespective of introducing 
the magnitude cutoff. 
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4. APPLICATION TO THE EXISTING LENS SURVEYS 
4.1. QSO sample and selection functions 

We apply the current lens model to the sample of 867 unduplicated QSOs at zs > 1 (as 
in Kochanek 1996), which are taken from the optical lens surveys such as the CFHT Survey 
(Crampton et al. 1992; Yee et al. 1993), the ESO Key-Programme Survey (Surdej et al. 1993), the 
HST Snapshot Survey (Maoz et al. 1993), the NOT Survey (Jaunsen et al. 1995), and the FKS 
Survey (Kochanek et al. 1995). The 5-band absolute magnitudes of these QSOs are evaluated 
from the V-band apparent magnitudes tabulated in each survey, a power-law QSO spectrum of 
oc z/~ 0 - 5 for the K-corrections, and a typical B — V = 0.2 mag in the QSO rest frame as adopted 
by previous workers. We note that V magnitudes of some numbers of QSOs are based on the 
catalogue of Veron-Cetty and Veron (1996) which provides not always accurate magnitudes. In 
this combined sample, only six QSOs are lensed (1208+1011, 1413+117, 1009—0252, 1115+080, 
0142—100, and 0957+561). We do not use 0957+561 having a large image separation Ad = 6."1, 
because the lensed image of this QSO is affected by an intervening cluster of galaxies. 

When we apply the lens model to these lens surveys, the selection functions S (see eq.3) 
inherent in the surveys for detecting multiple images give the additional constraints on the range 
of the impact parameter for lensing and thus affect the estimate of lensing probabilities. For 
each QSO, we evaluate the selection function which depends on the seeing in the relevant survey 
and the magnitude difference between the primary and secondary images. When the QSOs are 
observed in both the HST and other grand-based surveys, we adopt the HST results because of 
the high resolution of detecting separate images down to A 9 = ()" . 1 . 

In addition to the optical sample, we consider the radio lenses, as in Kochanek (1996). While 
the incomplete information on the redshift and 5-band luminosities of the sample in the radio 
surveys prevents us from calculating the absolute lensing probability in this sample, it is possible 
to use the observed image separations of the radio lenses in order to obtain how likely these data 
are reproduced in the current lens models. We adopt the ten radio lenses listed in Kochanek 
(1996) (CLASS 1608+656, MG 1131+045, MG 0414+0534, MG 1654+1346, B 1938+666, MG 
1549+3047, CLASS 1600+434, B 1422+231, MG 0751+2716, and B 0218+357). For the purpose 
of comparing his results with ours, we follow his procedure to assign redshifts and magnitudes for 
the lenses when these data are not available. 

We then calculate the sum of the lensing probabilities pi for the optical QSOs to 
obtain the expected number of lensed QSOs n = Jf,Pi and the image-separation distribution 
n(A9) = Jfpi(A9). The former is restricted to the lenses with A 8 < A" in order to avoid 
ambiguities associated with unusually large separations. We also calculate the likelihood function 
L for the models to reproduce Njj unlensed optical QSOs, Nr lensed optical QSOs, and Nr radio 
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lenses having image separations p{A0) (see Kochanek 1996), 

lni = -f> + |>*<Aty+|>(5^) , (11) 

i =1 j =1 fc=l \ Pk / 

where only the relative likelihood is computed for radio lenses to achieve the observed image 
separations. 


4.2. Results 

Figure 6 (a) shows the predicted number of optical lenses n(A9 < A") against Do for the flat 
universe of Do + Ao = 1. We take a standard parameter set of r* ore = 0.1/r -1 kpc, 7 = 4, rj = 1.2, 
and WN’s LF for QSOs. Different curves correspond to the results from the LFs of EEP (solid 
line), LPEM (dotted line), MGHC (dashed line), and K96 (thick solid line). For EEP, LPEM, and 
MGHC, three cases of M cut = —16.5, —17.0, and —17.5 mag are shown from the lower to upper 
curves, whereas the magnitude cutoff is not used for K96. Comparing the predicted number of 
lenses with the observed five lenses among the current QSO sample, it is obvious that different 
LFs predict the different best values of Do, as Do ~ 0.2 for MGHC, 0.3 for LPEM, and 0.5 for 
EEP. This effect of adopting different LFs overwhelms the effect of changing M cut significantly. 

It is remarkable that use of the updated LFs for E/SO galaxies yields a small Do (or large Ao) in 
contrast to the earlier strong conclusion in favor of Do = 1 (MR; K96). 

In order to elucidate the most likely case, we show in Fig. 6 (b) the image-separation 
distribution n(A9) for each LF with the value of Do that reproduces the observed number of 
lenses. The histogram denotes the image-separation distribution of the five optical lenses identified 
in the current sample, and the distribution of the ten radio lenses is shown by marking an asterisk 
at each image separation. Comparing the predicted and observed distribution n(A9) including 
radio lenses, the model using LPEM’s LF appears to reproduce the data better than the models 
using other LFs. Specifically, MGHC’s LF yields too many small separations and at the same time 
falls short of the observed number of both optical and radio lenses at 2" < AO < 3". This large 
difference between the predicted and observed separation distributions holds also for the LFs of 
both EEP and K96. 

Following these results, we determine the best combination of LF and cosmological parameters 
that maximizes the likelihood L defined in eq. (11) to simultaneously reproduce the number of 
lenses n(A8 < A") and the separation distribution n(A0). Figure 7 shows the likelihood results 
when using the LFs of LPEM (dotted lines), EEP (solid lines), and K96 (thick solid line), where 
three lines for LPEM and EEP correspond to the cases with M cut = —16.5, —17, and —17.5 mag. 
All likelihood values are normalized by its maximum L max which is derived using LPEM’s LF with 
M cut = —17 mag. Note that the case for MGHC’s LF yields quite a small L (In L/L max < —19) 
and is not shown in the figure. We have confirmed Kochanek (1996)’s result that an Do — 1 
universe gives the highest likelihood (In L/L max ~ —2.8) as long as his LF is used. However, EEP’s 
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LF with M cut ~ —17 mag also provides the similar likelihood at a lower Hq. It is evident from 
these cases that LPEM’s LF yields much larger likelihood over a broad range of Do- The peak of 
the likelihood is located at Do — 0.3 in flat cosmologies, and this appears to be very insensitive to 
the value of M cut . If we take Do = 0.3 as the most likely case, a universe with Do = 1 is ruled out 
at 98 % confidence level. For an open universe with Ao = 0, the likelihood monotonically increases 
with decreasing Do and In L/L max is as small as —1.5 at Do = 0.1 with LPEM’s LF. 


4.3. Uncertainties 

The above analyses support that the observations are best explained by the lens model using 
the updated E/SO-type LF of LPEM, among the LFs adopted in the current work. It should 
be kept in mind that MGHC’s LF contains some uncertainties related to the Zwicky magnitude 
system used in the CfA survey. We then examine how the likelihood based on LPEM’s LF is 
changed when the inherent model parameters are varied, and the result is tabulated in Table 3. 
Decreasing a core radius does not affect the likely range of Do significantly, because the effect 
of magnification bias is simultaneously decreased, as discussed in §2.2. Artificial increase of cr* 
irrespective of the value derived from the Faber-Jackson relation (eq. 9) tends to favor a larger Do, 
but the likelihood value turns out to decrease at the same time, as discussed below. Intriguingly, 
the likelihood result for LPEM’s LF is essentially unchanged even if no cutoff is employed, because 
of the shallow faint-end slope a. Use of MR’s QSO LF slightly decreases the likely value of Do, 
but the change is modest. 

The lensing probability is proportional to cr* 4 , thereby the statistics based on the observed 
number of lensed QSOs is sensitive to what value is assigned to cr*. We have estimated cr* from 
the Faber-Jackson relation (eq. 9) using the observed value of M* and the possible fraction of 
E and SO types. However, while the large fraction of SO-type galaxies for a given fraction of 
both E and SO types implies smaller cr*, as suggested from eq.(9), there are large uncertainties in 
distinguishing between E and SO types. Also, the self-consistent modeling of internal dynamics 
in E/SO galaxies suggests that cr increases slowly with increasing r core (Kochanek 1996). In view 
of these uncertainties in assigning the value of cr*, we in turn compute the likelihood for various 
values of a* in flat cosmologies, using LPEM’s LF. 

Figure 8 shows the likelihood in the two dimensional parameter space (cr*, Do) for flat 
cosmologies. Contours are shown at 68 % (1 cr), 90 %, 95.4 % (2 cr), and 99 % confidence levels 
for one degree of freedom in the likelihood ratio. It follows that within the 1 cr confidence interval 
for cr*, which includes cr* ~ 205 km s _1 derived from the Faber-Jackson relation (see table 1), a 
low-density universe with Do ~ 0.3 is favored. If we further increase a*, the peak of the likelihood 
is located at a higher Do but the overall likelihood turns to decrease. This is because the models 
with such large values of cr* yield larger average image separations than observed. 

Therefore, even if uncertainties in the lens model are taken into account, a low-density, flat 
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universe with a large Ao is statistically consistent with the observations when we use the updated 
LF for E/SO galaxies. In sharp contrast to the previous models of lensing statistics that have 
supported a high-density universe with Qq = 1, we conclude that a universe with Qq = 0.3^q'^ and 
+ Ao = 1 casts the best case to explain the results of the observed lens surveys. 


5. DISCUSSION AND CONCLUSIONS 

We have presented new calculations of gravitational lensing statistics in view of the recently 
revised knowledge of the internal dynamics and number density of early-type galaxies. Main 
revised points in the lens models are summarized as follows. (1) The factor of (3/2 ) 1 / 2 correction 
for a is not used following the detailed dynamical modeling of early-type galaxies (Breimer &; 
Sanders 1993; Franx 1993; Kochanek 1994). (2) The type-specific Schechter parameters for E/SO 
galaxies are adopted from the recent redshift surveys of galaxies combined with morphological 
classification. (3) The faint part of the LF of E/SO galaxies, which is dominated by diffuse 
dwarf-type galaxies, is found to be unimportant for lensing statistics. All of these revisions point to 
smaller lensing probabilities than earlier predicted. Applications of the lens models to the existing 
lens surveys suggest that both the total number of optical lenses and the image separations of 
optical and radio lenses are best reproduced in a low-density, flat universe. As the best set of 
cosmological parameters we arrive at S2o = 0.3 //} 2 with Hq + Ao = 1. A flat universe with Oo = 1 
is ruled out at 98 % confidence. We also find that Ao < 0.9 at 94 % confidence in flat cosmologies, 
and that an open universe with Ao = 0 is less likely than a flat universe with Ao 7 ^ 0 . 

Our conclusion supporting a large Ao is virtually in contrast to that of Kochanek (1996) who 
strongly argued against it using the same QSO sample as adopted here. Main difference between 
his result and ours originates from the points (2) and (3) stated above. In particular, Kochanek 
gathered the Schechter parameters of E/SO galaxies from various LF determinations by EEP, 
MGHC, and LPEM. He has adopted M* = —19.9 mag from the all-type LF by EEP, a = — 1 
from the all-type LF by LPEM and MGHC, and (/)* = 6.1 x 10 -3 from the scaling of LPEM’s 
(/)* = 14.0 x 10 3 for all types by MGHC’s E/S0 fraction. However, since the Schechter parameters 
have been determined by means of a highly correlated fitting in each galaxy survey, a simple 
mixture of these parameters taken from various references lacks consistency. Specifically, his use 
of brighter M* and steeper a than those of E/S0 galaxies reported in the original references (see 
Table 1) leads to a bias in favor of a systematically large lensing probability. A further effect 
against a large Ao seen in his model as well as other earlier models lies in the counting of large 
numbers of less compact dwarf-type galaxies as lens objects which dominate the faint part of the 
LF. 


Models for gravitational lensing must explain not only the observed probability of lensing but 
also the relative probability of showing a specified image separation. The average image separation 
increases with increasing a*, thereby the lens data are expected to constrain a* irrespective of the 
value derived from M* found in each galaxy survey and the Faber-Jackson relation. However, 
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our exercise of using the trivial effect of faint galaxies in lensing, which is also consistent with 
shallower faint-end slope a for earlier types seen in LPEM and LCRS (Lin et al. 1996; Bromley 
et al. 1997), implies a larger average image separation even when a* is fixed. As a result, the 
estimate of a* which accords with the image-separation distributions depends on the faint part of 
the LF for E/SO galaxies adopted. In this respect, it is notable that the lens model using LPEM’s 
a = 0.2 which was not explored by Kochanek (1996) provides much better fitting to the observed 
image separations for optical and radio lenses, even if his preferred value of a* — 220 km s -1 is 
used in our analysis. 

Intriguingly, our 1 c-confidence limit 0.2 < flo < 0.5 with Qq + Ao = 1 is in good accord 
with the recent result of Im, Griffiths, & Ratnatunga (1997) who analyzed the redshifts of known 
lensed QSOs combined with the observed image separations and lens properties. It is worthwhile 
to note that their method is immune to the uncertainties associated with the determinations of 
the LF. This may suggest that lensing statistics using the realistic form for the LF, as adopted in 
the present work, yield the similar results. However, before concluding so definitely, we require 
a sample of more high -2 QSOs in light of gravitational lensing, together with more definite 
knowledge of the internal dynamics and number density of lensing galaxies over a wide range of 
redshifts. Significant increase and improvement of the information on lensing phenomena will 
then allow us to answer whether there is a concordance of a flat, Ao-dominated universe with 
observations. 


The authors are grateful to T. Futamase and B. A. Peterson for useful discussions during 
the course of the present work. This work has been supported in part by the Grand-in-Aid for 
Scientific Research (09640328) and COE Research (07CE2002) of the Ministry of Education, 
Science, and Culture in Japan. 


A. The solutions to the lens equation 

From the lens equation (2), the critical impact parameter Iq for multiple imaging is written as 


IS = rl\(l + 5/3 - i/3 2 ) - + 4 ) 3 / 2 ] , (Al) 

provided (3 = r core /r e < 1/2. After some algebra, the three solutions to eq.(2) are given by 

h = -^ + 2p 1/2 cos^ (A2) 

b2 = -^ + 2p 1/2 cos(^ + y) (A3) 

b3 = ~\ l + 2 P 1/2 cos (^ - y) i ( A4 ) 
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Figure 9(a) shows the total amplification factor A = A\ + A 2 + A 3 for two representative cases 
of (3 = 0.1 and 0.25, compared to the SIS case. It follows that the presence of a finite core leads to 
diverging A at the caustics l = 0 and Iq, and hence results in the value of A which is larger than 
the SIS case. Figure 9(b) shows the cross section a = ttIq for multiple imaging (thick solid line), 
and the separation angles A 6 between two outer images (image 1 and 2) for l = 0 (solid line), 
0.7lo/r e (dotted line), and 0.9lo/r e (dashed line), respectively. It is clear that the presence of a 
finite core reduces the cross section for lensing, and that A 8 is almost independent of the impact 
parameter l, thereby being well approximated as y/\ — 2 (3 at l = 0. 
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TABLE 1 

E/SO Luminosity Functions' 1 and Lens Parameters 


Survey M* a (f>* a* F* 

(B mag) (Mpc -3 ) (km s^ 1 ) 

EEP 6 . -19.90 -1.10 4.8 x 10“ 3 214.5 0.018 

LPEM C . -19.71 +0.20 3.2 x 10“ 3 205.3 0.010 

MGHC d . -18.99 -0.90 9.1 x 10“ 3 166.0 0.011 

MGHC2 6 . -19.49 -0.90 9.1 x 10“ 3 186.2 0.017 

LCRS-f. -19.12 -0.27 1.1 x 10“ 2 179.2 0.017 

K96. -19.90 -1.00 6.1 x 10“ 3 225.0 0.026 


a In the B band and h = 1. 

b (j)* for all galaxy types is scaled by the E/SO fraction (31%) from 
Postman & Geller (1984). 

c 4>* is taken from the fit by Driver, Windhorst, & Griffiths (1995). 
a* is calculated from eq.(9) weighted by the Postman & Geller fraction of 
E (12%) and SO (19%) galaxies. 

d M* is an average of — 19.23B mag (E) and — 18.74B mag (SO), and 
a is an average of —0.85 (E) and —0.94 (SO), a* is calculated from eq.(9) 
weighted by MGHC’s fractions of E and SO galaxies. 

e M* is transformed from that of MGHC by 0.5 mag. 

/ For galaxies with no [O II] A3727 emission lines and with the mean color 
< bj — R >o= 1.1 mag. u* is from eq. (9) with the Postman & Geller fraction. 

References. -EEP, Efstathiou et al. (1988); LPEM, Loveday et al. 

(1992); MGHC, MGHC2, Marzke et al. (1994); LCRS, Lin et al. (1996); 

K96, Kochanek (1996) 
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TABLE 2 

Model Parameters for LFs of E/SO and dE/dSO Types' 1 * 


Galaxy Type M* a 4>* M cut 

(B mag) (Mpc -3 ) (B mag) 

E/SO 6 . -19.7 ^L0 3.2 x 10~ 3 -17.0 

dE/dSO. -16.7 -1.5 3.2 x 10” 3 

Both Types c . —20.0 —1.4 1.6 x 10~ 3 


a In the B band and h = 1. 

6 The Schechter function is multiplied by a factor exp(— io°- 4 * ( M_MCU< )). 
c Fit to a single Schechter form for the combined LFs of both E/SO 
and dE/dSO types. 
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TABLE 3 

Results of Maximum Likelihood for Optical and Radio Lenses 


Model 

Maximum 6 

68%(l(i) confidence 

90% confidence 

Standard Case a . 

. O 0 ~ 0.30 

0.17 < O 0 < 0.49 

0.12 < fi 0 < 0.66 

Parameter Change 

* 

r 

' core 

->• 0.05/r -1 kpc. 

. O 0 ~ 0.30 

0.19 < Q 0 < 0.54 

0.14 < fi 0 < 0.76 

r * 

' core 

->• 0.03/r -1 kpc. 

. ~ 0.35 

0.21 <n 0 < 0.59 

0.15 < fi 0 < 0.84 

a* -> 

210 km s -1 . 

. fi 0 - 0.30 

0.19 < n 0 < 0.53 

0.14 < Q 0 < 0.74 

a* -»• 

220 km s -1 . 

. O 0 ~ 0.40 

0.25 <n 0 < 0.65 

0.20 < fi 0 < 0.91 

M cut 

—* +oo mag. 

n 0 - o.3o 

0.18 < O 0 < 0.50 

0.13 < fi 0 < 0.70 

4>q ->• 

MR. 

fl 0 ^ 0.25 

0.13 < O 0 < 0.40 

0.05 <Q 0 < 0.56 


a Our standard choice of the parameters includes no (3/2) 1 / 2 correction for a*, 
LPEM’s <fig, WN’s 4>q, r* ore = 0.1 h~ l kpc, and M cut = —17 B mag, for a flat universe 
with Qq + Ao = 1. 

6 The value of Qq at L = L max when the binning interval for flo is 0.05. 
References. -LPEM, Loveday et al. (1992); WN, Wallington & Narayan (1993); 
MR, Maoz & Rix (1993) 
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M q (B mag) 

Fig. 1.— Magnification bias B defined in eq.( 8 ) as a function of QSO magnitudes Mq for h = 0.5. 
We adopt 7 = 4, 77 = 1.2, and a* = 220 km s -1 for the lens parameters, (a) Comparison between 
the SIS case (solid lines) and the case with a finite core r* ore = O.l/i -1 kpc for a = — 1 (dotted 
lines) and a = 0 (dashed lines). WN’s QSO LF is utilized. Note that presence of a core leads to 
large magnification bias and that this effect becomes larger for shallower faint-end slopes a of the 
LFs of galaxies, (b) Effects of using different QSO LFs on bias (WN: solid lines, MR: dotted lines, 
FT: dashed lines) for the case of r* ore = 0.1 h -1 kpc and a = —1. WN’s QSO LF appears to give 
the largest bias. 
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Fig. 2.— Lensing probability of a QSO with Mq = —27 mag (for h = 0.5) and zs = 2, as a 
function of image separations Ad. Various Schechter-form LFs of E/SO galaxies, as tabulated in 
Table 1, are examined, for a universe with (Qq,Xq) = (1,0). 
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Fig. 3.— Model LFs of early-type galaxies as a function of galaxy magnitudes M (for h = 1) 
are shown using the parameters given in Table 2. Dashed and dotted lines denote the LFs of 
E/SO and dE/dSO galaxies, respectively, and solid line denotes the combined LF of these galaxies. 
Approximate fit of the combined LF to a single Schechter form is shown by thick solid line. 





- 25 - 




q I_i_i_i_i_I_i_i_i_i_I_i tt i 1 ■ 

0 12 3 


Ad (arcsec) 

Fig. 4.— (a) Lensing probability at zl = 0.5 per unit galaxy magnitude dp/dM (normalized by 
F* for E/SO galaxies), when we use the model LFs shown in Fig.3. Calculations are made for a 
QSO with Mq = —27 mag (for h = 0.5) and zs = 2, and for (fio,-V)) = (1,0). See the text for 
the adopted lens parameters relevant to each galaxy type, (b) Image-separation distribution p{A6) 
predicted from the model LFs. The case for dE/dSO galaxies is not shown because of its giving 
negligible probability. 





-26- 



A 9 (arcsec) 



M cut (mag) 

Fig. 5.— Effects of introducing a cutoff magnitude M cut for the LFs of EEP (solid lines), LPEM 
(dotted lines), and MGHC (dashed lines), on the lensing probability. Calculations are made for 
a QSO with Mq = —27 mag (for h = 0.5) and zs = 2, and for (flo^o) = (1,0). (a) Image- 
separation distribution p(A9) for M cut = +oo (no cutoff), —16, and —17 mag (for h = 1), and for 
r core = O.l/i- 1 kpc. (b) Ratio of the total lensing probability p with cutoff relative to that without 
cutoff, as a function of M cut . Two cases for r* ore = 0.05 and 0.1 h^ 1 kpc are shown. 
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Fig. 6.— (a) Predicted total number of lenses n with A 9 < A" in the adopted optical lens surveys, 
compared with the observed five lenses (thin solid line). Model calculations are made for a flat 
universe Do + Ao = 1, and other parameters are given in the text. Different model curves correspond 
to the results from the LFs of EEP (solid lines), LPEM (dotted lines), and MGHC (dashed lines), 
with cutoff magnitudes M cut = —16.5 (upper lines), —17 (middle lines), and —17.5 mag (lower 
lines), respectively. For comparison, the case for K96’s LF without employing cutoff (thick solid 
line) is also shown, (b) Predicted image-separation distribution n(A0), compared with the observed 
image-separation distribution in the optical sample (histogram) and in the optical lenses (asterisks 
located at their respective separations A 9). The value of Do for each model with cutoff M cut = —17 
mag (except for K96) is chosen so as to reproduce approximately the observed number of optical 
lenses, examined from panel (a). 
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Fig. 7.— Results of the maximum likelihood analyses for reproducing both the total number of 
optical lenses n with A 6 < 4" and the image-separation distribution n(A9) of optical and radio 
lenses. The likelihood functions L defined in eq.(ll) are shown as a function of Do for a flat 
universe Do + Ao = 1. Different curves correspond to the results from the LFs of EEP (solid 
lines), LPEM (dotted lines), and K96 (thick solid line). For EEP’s and LPEM’s LFs, three cases 
of M cut = —16.5, —17, and —17.5 mag are shown by lower, middle, and upper lines at the value of 
Do = 0.1, respectively. All likelihood values are normalized by its maximum L max which is derived 
using LPEM’s LF with M cut = —17 mag. The confidence levels are indicated on the right-hand 
side of the plot. Note that LPEM’s LF at Do — 0.3 yields by far the large likelihood compared 
with other cases. 




Fig. 8.— Likelihood contour plots for flat cosmologies in the two dimensional parameter space 
(<t*,S7o)) f° r our standard model using LPEM’s LF with M cut = —17 mag. Contours are shown 
at 68 % (1 a), 90 %, 95.4 % (2 a), and 99 % confidence levels for one degree of freedom in the 
likelihood ratio. 




- 30 



Fig. 9.— (a) Total amplification factor A = A\ + A 2 + A 3 of three lensed images as a function of 
the source position l normalized by r e . Solid lines correspond to the cases for non-zero core radii, 
(3 = r core /r e = 0.25 and 0.1, while dotted line for the SIS case shows the amplification factor of two 
lensed images, (b) Ratio of the lensing cross section a = ttIq to that for the SIS asis = Trr’e (thick 
solid line) and the ratio of the image separation A 6 to that for the SIS A Osis (solid, dotted, and 
dashed lines), as a function of f). 





